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We study real space condensation in aggregation-fragmentation models where the total mass is 
not conserved, as in phenomena like cloud formation and intracellular trafficking. We study the 
scaling properties of the system with influx and outflux of mass at the boundaries using numerical 
simulations, supplemented by analytical results in the absence of fragmentation. The system is found 
to undergo a phase transition to an unusual condensate phase, characterized by strong intermittency 
and giant fluctuations of the total mass. A related phase transition also occurs for biased movement 
of large masses, but with some crucial differences which we highlight. 
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Condensation transitions constitute an impor- 
tant class of non-equilibrium phase transitions, and 
occur generically in many mass transport models Jl[ 
such as the zero range process and its variants 
and the aggregation-chipping model These sys- 
tems are characterized by a fixed total mass (num- 
ber of particles) and stochastic rules for exchange of 
mass between sites. When the total mass of the sys- 
tem exceeds a critical value, condensation sets in, 
with a finite fraction of the total mass forming a 
macroscopic cluster that occupies a single site. The 
phenomenon is akin to Bose condensation, but in 
real space. 

Does the condensation transition survive in a 
system when the total mass is not conserved, but can 
undergo large fluctuations due to the exit of clusters 
of all sizes? This question is important in a num- 
ber of physical situations, ranging from formation of 
clouds and aerosols, to intracellular trafficking and 
organelle formation in living cells 0-01 • We address 
this within a simple but generic ID model with ag- 
gregation and fragmentation (chipping) of masses in 
the bulk, and influx and outflux of masses at the 
boundaries. Our main flnding is that the open sys- 
tem does undergo a condensation transition upon in- 
creasing the influx or decreasing the chipping rate. 
However, the nature of the condensate is very dif- 
ferent from that in the closed model Q, in that the 
mass in the condensate shows giant number fluctu- 
ations and has a broad distribution, in contrast to 
the sharply peaked distribution in the closed system 
0, 0i- The condensate, however has a well-deflned, 
finite mean mass for a fixed system size and is thus 
quite different from the indefinitely growing aggre- 



gates in open models which allow only single parti- 
cles to exit at the boundaries Q . 

The intermittent and fiuctuating nature of the 
condensate gives rise to novel signatures: the to- 
tal mass M itself shows giant fiuctuations and has 
a distribution characterized by a prominent non- 
Gaussian condensate tail whose width scales with 
system size. Further, the exit of the condensate from 
the boundaries and the accompanying sharp drops 
in M give rise to interesting 'charge and fire' be- 
haviour of M: the time series M{t) departs strongly 
from self-similarity and shows quantitative features 
of intermittency, which we characterize in terms of 
appropriately defined structure functions, as in tur- 
bulence phenomena. Turbulence, in the sense of 
multi-scaling of n-point mass-mass correlation func- 
tions has been studied earlier in aggregation mod- 
els 0, but our characterization of turbulence-like 
behaviour is quite different, being associated with 
temporal fluctuations of total mass. Our results are 
based on both analytical and numerical work. In the 
limit of zero chipping, we analytically calculate the 
moments of total mass in steady state, and also the 
dynamical structure functions, whereas for non-zero 
chipping, we perform numerical Monte Carlo simu- 
lations. 

Recently, it has been demonstrated that giant 
number fluctuations are related to anomalous, non- 
Porod behaviour of spatial correlation functions in 
a wide class of systems [l^. Our work points to a 
connection between giant number fluctuations and 
anomalous dynamical behaviour, namely temporal 
intermittency, which is related to higher order cor- 
relation functions in time 1 121. It also raises the 



interesting general question of whether temporal in- 
termittency is present in other systems with giant 
number fluctuations and suggests that dynamical 
structure functions, as used in the paper, provide a 
useful probe of intermittency in these systems also. 

We work with a general lattice model incorpo- 
rating diffusion, aggregation, fragmentation, influx 
and outflux, which goes beyond earlier studies of ag- 
gregation with input |13l4l5l | , and aggregation and 
fragmentation in a closed system [3|, |l6|. Starting 
with an empty lattice of L sites at t = 0, a site i is 
chosen at random, and one of the following moves 
occurs: 



forward and backward chipping rates are unequal, 
an aggregate is not expected to form [l7|. Thus, 




FIG. 1. Model: Influx of unit masses at site 1. Forward 
and backward stack movement at rates D and D' respec- 
tively. Forward and backward chipping at equal rates w. 
Outflux of full stack or unit mass (via chipping) from 
site L (and 1). 



i. Influx: A single particle of unit mass is injected 
at rate a at the first site {i = 1). 

ii. Diffusion and aggregation: With rate D (or 
£)'), the full stack on site i (i.e. all particles 
on the site collectively) hops to site i + 1 (or 
i — 1) and adds to the mass already there. 

iii. Chipping of unit mass: With rate 2w, a unit 
mass breaks off from the mass at i and hops to 
site i — 1 or z + 1 with equal probability, adding 
to the mass already there. 

iv. Outfiux of mass from boundaries: With rate 
D (or D'), the entire mass at site L (or site 
1) exits the system; with rate w, a unit mass 
breaks off from site L (or site 1) and exits. 

We find that the results depend strongly on two 
factors: one, whether motion of particles is biased or 
not and two, whether or not exit of masses is allowed 
from the boundary where influx occurs. In this pa- 
per, we only consider the effect of bias IT]. We find 
that the occurrence of the phase transition is robust 
with respect to bias in the movement of stacks, but 
not chipping. As in the closed periodic case, if the 



chipping is taken to be unbiased in both the cases 
we study in this paper: 

(A) Unbiased Stack Hopping: D = D'; exit allowed 
from sites 1 and L. 

(B) Biased Stack Hopping: D' — 0; exit allowed 
from site L 

Infiux and chipping occur at rates a and 2w respec- 
tively in both the above. 

(A) Unbiased Stack Hopping (£»' = D): We dis- 
cuss both the phases and the critical point below: 

Normal (large w ) phase: In this phase, a typical 
configuration does not show very large fluctuations 
about the average mass profile. The total mass M 
too has normal fiuctuations i.e. the rms fluctuations 
AM = ,/(M2) - (M)2 cx a/I, with the distribu- 
tion for the rescaled mass variable {M — (M))/AM, 
approaching a Gaussian at large L. The mass distri- 
bution P(m, j, L) at a given site j is found to depend 
on j and L only through the rescaled position vari- 
able X = j/L |18| . implying that for a given x, all 
moments of mass are independent of L to leading 
order. 

Condensate ( small w ) phase: A typical config- 
uration deviates strongly from the average profile, 
with the largest (local) fiuctuations scaling as sys- 
tem size L. On monitoring the largest mass mi 
in the system, we find that its average value (mi) 
is proportional to L implying that the system 
contains a macroscopic condensate. The presence of 
the condensate has a strong effect on all steady state 
properties of the system such as the total mass M, 
mass at a site, etc. The probability distributions of 
all these quantities have an exponential tail with a 
characteristic mass Mg where Mq oc L for a given 
w and a. We refer to this exponential tail as the 
'condensate' tail and describe below, how it appears 
in various steady state distributions: 

i. The steady state distribution P{M) of total 
mass of M in the system behaves as P{M) ^ 

^exp(-^) at large M [fig. [2] . Conse- 

quently, the rms fluctuation AAf of total mass 
shows non-Gaussian behaviour, scaling as L 
rather than ^/L. We have analytically cal- 
culated various moments of the total mass in 
the limit w = [HI. We find that AM/L ~ 
0.46(a/D), in the limit of large i, which agrees 
well with numerics. 

ii. The distribution P{mi) of the largest mass mi 
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FIG. 2. P(M) vs. M for L = 100 and L = 200 in the 
normal phase {a — 1, D = 0.75, = 2) and conden- 
sate phase (a = 1, D = 0.75, w = 0.25). Inset: Scahng 
collapse of tails on plotting LP{M) vs. M/L in the con- 
densate phase and L^'''^P{M) vs. {M - {M))/L^''^ near 
the critical point {a = 1, D = 0.75, w — 1.5). 



also follows P(mi) ^ exp (^^^^ for large 
mi [isi]. 

ill. The distribution of masses exiting from the 
left or right boundary [3| is found to follow: 
Pexttim) - (ilo exp (-i^j), for large m 

iv. The single site mass distribution P{m,j,L) 

m follows i/ ii) exp (-^)) at large 

m [l^ . The factor 1 /L arises as the aggregate 
can be at any one of the L sites, and f{j/L) 
reflects that the aggregate does not visit all 
sites with the same probability. The rms fluc- 
tuation Am{x,L) of mass at a given x = j/L 
is thus anomalously large as well: it increases 
as y/L with L rather than being 0(1), as in the 
normal phase. 

That there is no constraint on the total num- 
ber of particles per site in our model is crucial for 
i-dependent fluctuations to arise. Systems such as 
vibrated needles [20*1 and passive particles in fluctu- 
ating fields [21i . .22] also display giant number fiuctu- 
ations, but in these systems, fluctuations in a region 
of linear size Al depend primarily on Al, rather than 
L [l3| • This is traceable to hard core interactions be- 
tween particles in these models. Once this constraint 
is removed, macroscopic stacks can form and mass 
fluctuations depend on L ■ 

Critical point Wc-' The transition from the nor- 
mal to the condensate phase takes place at a crit- 



ical chipping rate Wc, which increases with injec- 
tion rate a if D is held constant At criticality 
also, large fluctuations of the total mass are found, 
consistent with AM cx L^" where 9c — 2/3. The 
mass distribution has a tail of the form: P{M) ~ 

_i_ exp ( "^^ M^^^^O ^^^^^ ^^'^^ ^ ^ ^^'^ ^'^^ ^ 
[inset, flg. [2] . Interestingly, we flnd that there 
is a similar L-dependent tail in the distribution of 
masses exiting from the left, but not the right. This 
is presumably because although an L-dependent ag- 
gregate forms close to the left boundary, it dissipates 
due to chipping on diffusing through the bulk of the 
system, and does not survive up to the right bound- 
ary. 

Contrasting signatures of the phases also ap- 
pear in dynamical properties: M{t) is self-similar in 
time in the normal phase but exhibits breakdown of 
self-similarity in the condensate phase. The break- 
down of self-similarity is captured in the behaviour 
of the structure functions S'„(t) = ([M(t) - M(0)]") 
25| where (...) denotes average over histories. Self- 
similar signals typically show Sn{t) oc as t/r — > 
0, where 7 is a constant and t is a time scale which 
characterizes the lifetime of the largest structures in 
the system. A deviation from Sn{t) oc t"^" reflects 
the breakdown of self-similarity and may occur, for 
example, if the signal M{t) alternates between peri- 
ods of quiescence (small or no activity) and bursts 
(sudden large changes) [l2j. Such an alternation 
is characteristic of intermittency. The most well- 
studied measures of intermittency are the flatness, 
deflned as K{t) — S'4(t)/5'|(t); and the hyperflatness 
h{t) — S'6(t)/S'|(i). For intermittent signals, both 
K{t) and h{t) diverge as t/r — )• [13]. Below we 
present evidence for intermittency in our model. 

Normal Phase: In this phase, the structure func- 
tions are independent of L at small t and scale as 
S2n ~ t^'^ where the dependence of /3„ on n is close 
to linear indicating self-similarity of the time 
series M{t) [fig. 3(a)| . The flatness K{t) and hyper- 
flatness h{t) approach a finite, L independent value 
ast^O [fig. [3(c)]and fig. 4 in [T8|]. 

Condensate Phase: In the condensate phase, 
M{t) builds up as mass is injected and drops as 
masses exit, with occasional large crashes [fig. |3(b)| 
corresponding to the exit of condensates with 0{L) 
mass. The structure functions are found to scale as: 
Sn ^ L""- fn{t / L"^) , where /„ is consistent with the 
form fn{y) ^ (-l)"yg„[log(y)l for smaU y, (with g„ 
chosen to be a polynomial) [26|, and ap pro aches an 
n-dependent constant value at large y [l8|. Thus, 
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FIG. 3. (a)-(b): Realizations of M{t) vs. t for different 
L in the (a) normal phase {a = 1, D — 0.75, w — 3.0) 
and (b) condensate phase {a = 1, D = 0.75, w = 0.25). 
Note that the y-axis in (a) and (b) has a different scale, 
(c) K{t) vs. t for L = 100 and L = 200 in the two phases. 
Solid lines are fits to the form described in the text for 
t <^ for w < Wc- Inset: Scaling collapse of n{t) vs. t 
for different L on scaling time as t/L'^ in the condensate 
phase. 



the system shows strong intermittency: at small t, 
all structure functions 5„ behave as ~ t with the n- 
depcndence entering only through the multiplicative 
\ogt terms. It follows that K,{t) and h{t) diverge at 



small times in a strongly L dependent way [fig. 3(c) 
and fig. 4 in [18]]. In fact, they are functions of tjl7 
and diverge as t/L"^ [inset, : 
fig. 4 in |18|]. We have also analytical! 



3(c) and inset, 
calculated 



S2{t) in the zero chipping limit w = Q [18j and find 
that it agrees well with numerical results. 

Critical point Wc-' M(t) continues to show inter- 
mittency at the critical point with fiatness and hy- 
perflatness diverging as i in an L dependent 
manner. However, there seems to be no simple scal- 
ing which collapses the curves for different L. 

(B). Biased stack hopping {D' — 0): The steady 
state can be obtained exactly in the limiting cases of 
only chipping I? = @ and only aggregation w = 
27[. The probability distribution of the rescaled 



mass, {M - (Af))/AM, is Gaussian with AM oc 
\/L, in both limits but for different reasons. In the 
pure chipping limit D — 0, this follows from the 
independence of masses at different sites, implied by 
the product measure of the mass distribution [§] ; in 
the pure stack hopping limit w = 0, on the other 
hand, it is associated with the formation and exit of 
aggregates of typical size ^/L . This essential 

difference is well captured by the time series data. 
For w = 0, the total mass M shows intermittency on 
time scales of order \/L, corresponding to a typical 
time interval of 0{VL) between exit events. Flatness 
and hyperffatness are functions of t/^/L and diverge 
as power laws as t/^L — >■ 0. By contrast, for D = 
0, the time series M{t) is not intermittent. Thus, 
intermittency rather than anomalous steady state 
fluctuations of M, is a key signature of aggregate 
formation when stack hopping is driven. 
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FIG. 4. {m^{x)) vs. X ioT a = 1, D = 1.5 and different 
w. Note the upward (downward) bending of curves on 
log-log plot in aggregation (chipping) dominated phase. 
There is no bending at Wc- 

As w is decreased, there is a phase transition 
from a normal phase to an aggregation-dominated 
phase characterized by intermittency. Unlike the un- 
biased case, however, the typical size of aggregates 
that exit the system is now expected to scale as -s/Z 
rather than L. This is consistent with the behaviour 
of (to^(x)) vs. X [fig. 11]. For large w, the plots of 
{m?(^x)) approach a constant value at large a;, thus 
indicating that there are no x dependent aggregates 
at large x and no intermittency. For small w, the 
plots of {m?(x)) vs. X bend upwards, consistent with 
an approach to ^/x at large x. Exit of Vl sized ag- 
gregates leads to intermittency as for w = 0. 
The transition takes place at Wc, (corresponding to 
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the curve with no bending on the log-log plot), at 
which {m?{x)) behaves as ~ a;"'= with etc — 0.16. 
M{t) shows intermittency at the critical point also. 

In conclusion, the principal result of this work 
is to establish the existence of a condensate phase 
in unbiased aggregation-chipping models where total 
mass is not conserved due to influx and outflux at the 
boundaries. This phase is characterized by anoma- 
lous steady state fluctuations of the total mass, and 
by intermittency in the dynamics, as quantified by 
the divergence of the flatness. It is likely that flat- 
ness would be a useful measure in other mass ex- 
change models also. The phase transition also occurs 
when the movement of stacks is biased, but the inter- 
mittent, aggregation-dominated phase in this case is 
different. 

Acknowledgements: We thank D. Dhar for useful 
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Condensation and Intermittency in an Open Boundary 
Aggregation-Fragmentation Model— Supplementary Material 



In this supplement, we provide the foUowing details of analytical and numerical results for the 
model with unbiased stack hopping: (i) an outline of the calculation of (M") and S2{t) in the limit 
m = of this model (ii) further numerical evidence for the differences between the condensate and 
normal phases, as described in the main paper (iii) a theoretical estimate of the critical chipping 
rate Wc- We also provide numerical evidence for intermittency in the aggregation-dominated phase 
of the model with biased stack hopping. 



UNBIASED STACK MOVEMENT 
A. Some analytical results in the limit to = 0: 

The condensate phase extends in the range < 
w < Wc- In the zero chipping limit w — 0, where 
stack movement is the only dynamical move in the 
bulk, we are able to obtain some analytical results 
for the moments of total mass in the steady state, 
and the structure function 52 (t). These calculations 
are outlined below [1]. 

Calculation of moments (M") in steady state: 

To calculate moments of total mass M , we consider 

i 

mass between sites i + 1 and j , i.e. mij — ^ 

k=i+l 

(where is the mass on the fc*'' site) and the cor- 
responding probability distribution Pij{m) of this 
mass. In steady state, Pij{m) satisfies the following 
equation; 

- iPi j{m) = < z, i + 1 < j, j < L 

(la) 



a[l - S,n,a]Po,]{rn - 1) - aPQj{m) + D[Poj+i{m) 
j_i(m) + Pi.jim) - 3Po,, (m)] - 



(lb) 



P^X-l{m) + P.-i,L(m) + P.+ix(m) - 3P,,L(m) = 
< z, i + 1 < L 

(Ic) 

2<5m,o + ^'i,i+2(m)+P,_i,,+i(m) -4P,,i+i(m) = 
0<i, 1 + 1= j<L 

(Id) 

To obtain (M"), we now do the following: 
(i) obtain recurrence relations satisfied by the gen- 

oo 

crating function Qij{z) ~ ^ Pij{m)z"^ from the 

m— 1 



above equation (ii) differentiate Qij{z) n times to 
get recurrence relations for the n*'' mass moments 
(ni^j) (iii) convert the recurrence relations satisfied 
by ) to a differential equation satisfied by (m"y) 
by taking the continuum limit x — i/ L, y — j / L (iv) 
solve this differential equation (which turns out to 
be a Laplace equation on a triangular region) to get 

to get 



"xyl 



(v) substitute a; = 0, y = 1 in 



(M") . This gives the following expressions for (Af") , 
to leading order in L: 



(M) = 0.5 



aL 



AM = Vp5^(M)F) - 0.459 



(M^) - 0.461 , 
aL 



(M^) - 0.615 



aL 
~D 



(M^)- 1.074 



(2) 



Calculation of structure function S2 (t) : 



S2{t)^{[M{t)-M{0)r) 

= 2[(Af2) - (Af)2] - 2[{M{0)M{t)) - (Af)2] 



(3) 



We have already calculated {M^) and (Af). To 
obtain the auto-correlation function (Af(O)Af(t)) — 
(Af)^, we calculate the quantity Gij{t) — 
(Af(O)mij (i)) — {M){mij). This satisfies the equa- 
tion; 



dt 



--D[{G,+i,j{t) + G,_i.j(t) + G,j+i(i) 



-HG.,,_i(i)-4G..,(i))] 



Go,,(i) = Gi,,(i) 
G,,.(t) = 



G. 



L + l 



(4a) 



(4b) 



G,,,(0)-(Afm,,,)-(M)K,j) 



(M)(m,j) 



'■OA/ 



3,L/ 



(4c) 



As before, we take the continuum limit x = i/L, 
y = j / L m. space to get a differential equation which 



we solve to obtain Gxy(t) and hence 82(1). 
gives the following expression for S2{t): 



E 

n=l,3,5. 



16?;^ / 

7 TT 

(n7r)4 V 

X < 1 — exp 



coth 



2 



- 1 



This 



8v 



(5) 



To extract the small t behaviour, we convert the 
above sum into an integral and take the limit r = 
^ — > 0. In this limit, the coth term in the above 
sum behaves asymptotically as — Tlog[r]. Thus, we 
obtain the following small t expression for S2{t) 



. , , . Dt 



(6) 



where ^0 f^nd Ai are some constants. Higher order 
structure functions can also be calculated similarly, 
but the calculations become very cumbersome. 



B. Numerical evidence 

The condensate and normal phases in the 
model are characterized respectively by the pres- 
ence or absence of a macroscopic condensate with 
typical mass which is 0{L) for a system of L sites. 
Signatures of this condensate appear in both the 
steady state and dynamical properties of the system. 
Below, we briefly recapitulate these properties, as 
presented in the main paper, and provide additional 
numerical evidence for them. 



Steady state 

In the condensate phase, the following steady 
state mass distributions have an exponential tail 
with characteristic mass Mq oc L, which we refer 
to as the condensate tail: 
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1. Exit currents: In the condensate phase (a = 1, 
0.75, w = 0.25), Pl/jitim) show condensate tails 



with good scaling collapse to Pl/jui'm) vs. m/L for 



different L. Inset: In the normal phase {a — 1, D — 0.75, 
w — 3.0), exit currents do not have condensate tails. 
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FIG. 2. Single site mass distributions: In the condensate 
phase {a — 1, D — 0.75, w — 0.25), P{m,j,L) has con- 
densate tails with the same characteristic A/q for differ- 
ent j/L. There is good scahng collapse of P(m,j, L) 
vs. m/L for a given j/L. Inset: In the normal phase 
(a = 1, D = 0.75, w = 3.0), P{m,j,L) for a given j/L 
is independent of L to leading order and shows no con- 
densate tail. 



i. P{M): The probability distribution P{M) of 
the total mass M in the system obeys P{M) ^ 

-pj^exp^— for large M where Mq oc L. 
Evidence for this is shown in fig. 2 of main 
paper. 

ii. Pexit{m): The distributions Pixiti''^) ^"^^ 
Plxit{ni) of masses m exiting from the 

sites 



left and the right boundaries (i.e 
1 and L) respectively obey 
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FIG. 3. Probability distribution of largest mass mi in 
the condensate phase (a = \, D = 0.75, w — 0.25): 
P{mi) has condensate tails with good scaling collapse 
to LP{mi) vs. mi/L. Inset: (mi) scales linearly with 
L. 



(iBo ^''^P (^^)) large m where the 
prefactor Ai/^ is different for the two distribu- 
tions. Figure [1] shows the scahng collapse of 
the tails of the two distributions on plotting 

L^p1\M) vs. M/L. 

iii. P(m, j, L): The mass distribution at the 
j'*^ site of the system obeys P{'m,j,L) ~ 

r/(i)(4exp(-^)) at large m. Fig- 

ure [2] shows the scaling collapse on plotting 
L^P{m,j/L = 0.4) and L'^P{mJ/L = 0.8) vs. 
m/L. As implied by the form of P{m, j, L), 
for large m, the single site distributions for 
different j/L collapse on to different curves, 
which differ from each other only by the con- 
stant prefactor f{j/L). 

iv. P(mi): The mass distribution of the largest 

mass mi behaves as P{mi) ^ exp (^—-^^ 

for large mi [fig. |3]. As a result, the largest 
mass in the system is macroscopic (on an av- 
erage), scaling as (mi) oc L [inset, fig. [3] 

By contrast, in the normal phase, we have: 

i. P{M): The distribution for the rescaled mass 
variable {M — (A/)) /AM, approaches a Gaus- 
sian at large L (see fig. 2 of main paper). 

ii. Pexiti'm): To leading order in L, Pl^^^{m) is 
independent of L, while Pl^^^{m) ~ {h/L)"^ 
where 6 is a constant. (See inset in fig. [T]). 
This behaviour of the exit currents can be 



established analytically for the limiting case 
D = 0. 

iii. P{m, j, L): The mass distribution P{m,j,L) 
at a given site j depends on j and L only 
through the rescaled position variable x = 
j/L. The inset in fig. [2]shows P{m,j/L = 0.4) 
and P{m,j/L = 0.8) vs. m for different L. 
Note that there is no significant dependence 
on L for a given j/L. 

iv. (mi) is sub-linear in L, such that {mi) / L — !■ 
for large L. 



Dynamical properties 

As discussed in the main paper, a key signature 
of the condensate phase is intermittency in the time 
series of the total mass M . This is captured by the 
small t behaviour of the structure functions Sn(t): 

Condensate phase: Numerics show that the 
structure functions scale with L as: S'„ ~ 
L^^fnit/L^)- Motivated by the small t form for 
5*2 (t) in the limit w = we use the form 
fn ^ (-l)"y5„[log(?/)] for small y, and get good 
fits with gn chosen as polynomials. Figure 4(a) 
shows the scaling collapse of the n*^^ order struc- 
ture functions for two different L to the above form 
for n = 2,4,6. We plot S2n/L^ vs t/L^ (where the 
rescaling L ~ i/100 has been done to just display all 
three structure functions in the same plot clearly). 
The solid lines are fits to tg2n[^og{t/L'^)] where 172, 
34 and ge are taken to be polynomials of degrees 1, 
2 and 3 respectively. The above forms of S'„ im- 
ply that the flatness K{t) = 5'4(t)/5'|(t) and hyper- 
flatness h{t) = Se{t)/S2it) diverge at small t. The 
divergence of the flatness in the condensate phase 
is shown in fig. 3(c) of the main paper. Figure [5] 
shows the divergence of the hyperflatness, with the 
solid lines obtained from the fitted forms of 5*6 and 
5*2. The inset in fig. [S] shows the scaling collapse of 
h{t) for different L on scaling time as t/L'^. 

Normal phase: At small times (t <C i^), the 
structure functions Sn{t) are independent of L and 
behave as Sn ~ i^" (except very close to t = 0), 
with the dependence of /3„ on n being close to lin- 
ear. Figure [4(b)] shows 52, ^4 and 5*6 for two dif- 
ferent L along with solid lines t^" with f32 ~ 0.52, 
f34 ^ 1.03 and (3^ ^ 1.45. The small deviation from 
Pn (X n may be due to sub-leading correction terms. 
Thus, in this phase, flatness (see fig. 3(c) of main 
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FIG. 4. Structure functions: (a) In the condensate phase 
(a = 1, D = 0.75, w = 0.25), structure functions for 
different L show good scaling collapse to S2n{t)/L^ vs 
t/L^ . Solid lines indicate tg2n\log{t / L^)] as described in 
text. 

(b) In the normal phase (a = 1, D = 0.75, w = 3.0), 
S2n{t) show no dependence on L and behave as ~ t^^"- 
(solid lines) for small t. 




FIG. 5. h{t) vs. f for L = 100 and L = 200 in the 
condensate (a = 1, _D = 0.75, w — 0.25) and normal 
{a = 1, D — 0.75, w = 3.0) phases. Solid lines are fits 
to the form described in the text for t <^ for w < Wc- 
Inset: Scaling collapse of h(t) vs. t for different L on 
scaling time as t/L^ in the condensate phase. 



We obtain a rough analytical estimate for Wc as fol- 
lows: the system is in the condensate phase as long 
as there is at least a finite fraction of the lattice that 
locally satisfies conditions for aggregate formation. 
As in the closed periodic case [l] , we take these con- 
ditions to be: (i) the occupation probability of a site 
in this region is equal to a critical occupation proba- 
bility Sc(w;) and (ii) the mass density is greater than 
a critical mass density Pc{w). While we do not have 
the forms of Sc and pc in the open boundary case, 
in the limit L — > 00, we approximate these by the 
exact conditions derived for the closed periodic case 



Sc{w) ^ iVTT^ -l)/iVTT^+l) (8a) 



paper) and hyperflatness (fig. [5]) approach finite, L 
independent values as i — >■ 0. 



C. Estimation of critical point Wc 

Let us define pj and Sj respectively as the aver- 
age density and occupation probability (probability 
that the site is non-empty) at a site j. In steady 
state, the balance of average mass current at each 
site gives the following exact expression (with D set 
equal to 1): 

pj+wsj =a{l-j/L) (7) 



Pj > Pc{w) where Pc{w) — yl + w — 1 (8b) 

Numerical simulations suggest that this is a reason- 
able approximation, at least away from the right 
boundary x — L. At Wc, a vanishingly small frac- 
tion {j /L — 0) of the system satisfies conditions 
for aggregate formation. Thus, setting Sj = Sc{w), 
Pj — Pc{w) and j /L in ([7]), we get: 

= (2a - 1 + \/4a + l) /2 (9) 

This agrees with the numerically determined Wc, ob- 
tained from AA/ vs. L curves, to within 10%. 
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tcrmittency of M{t), as the steady state mass dis- 
tribution P{M) has the same (Gaussian) form in 
both phases, as explained in the main paper. Fig- 
ured shows the flatness K{t) for two different values 
of L in both the aggregation-dominated and nor- 
mal phases. In the aggregation-dominated phase, 
the flatness K{t) diverges as i — in an i de- 
pendent manner while it approaches a constant L- 
independent value in the normal phase. The inset in 
fig. [5] shows the scaling collapse of K{t) for different 
L on rescaling time as t/^/L. Thus, intermittency 
occurs on time scales of order \/L when stack hop- 
ping is biased. 



FIG. 6. Biased stack Movement: K{t) vs. t iox L — 200 
and L — 400 in the aggregation-dominated {a — 1, 
D = 1.5, w = 0.125) and normal (a = 1, D = 1.5, 
w — 3.0) phases. Inset: Scaling collapse of it{t) vs. t for 
different L on scaling time as t/y/L in the aggregation- 
dominated phase. The divergence at small t roughly fol- 
lows ~ (t/\/T,)~^ (solid line). 

BIASED STACK MOVEMENT 

When stack movement is biased, the basic signa- 
ture of the aggregation-dominated phase is the in- 
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